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ABSTRACT 



This investigation is an extension of the work of Wrigley 
and Hovorka in fire control to the problem of the thesis. Quanti- 
ties peculiar to the underwater problem are defined, examined, and 
illustrated when appropriate. Using these quantities and those 
from the "above-water" problem that apply, the problem is formu- 
lated. Two possible methods of solution are discussed and 
compared . 
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CHAPTER I 



INTRODUCTION 



1 . 1 The Advent of the Stand-Off ASW Weapon 

A common conception of surface ship antisubmarine warfare 
is the picture of a destroyer plunging through the seas, dropping 
depth charges and firing short-ranged hedgehogs at a somewhat 
elusive target. This is not true today and will be even less true 
in the future. According to Vice Admiral John W. Thach, USN , 
Commander Antisubmarine Warfare Force, U.S. Pacific Fleet in a 
recent article; 

The new ASW destroyers will be very different from 
their predecessors. Their detection capability will 
be much longer ranged and sharper than it is today... 
Multi-purpose missile systems with accurate fire con- 
trol to provide for early attack at long ranges will be 
forthcoming in the next decade...' ' 

Thus the surface- to-underwater fire control problem is one of 

interest and its scope is changing rapidly. In the past, ASW 

ships little more than threw explosive rocks into the water when 

they thought they were about over the target. Even the throwing 

was slow and the time that elapsed between the order to throw and 

the throwing had to be taken into consideration. Such will not be 

the case in the future. Detection systems and weapons will permit 

a destroyer to engage his foe at long range. 

^'^The numbers in parentheses refer to cited references in 
the Bibl iography . 
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1.2 Previous Work in Fire Control 

This thesis is based on the work done in fire control 

theory at the Massachusetts Institute of Technology by Walter Wrigley 

( 2 ) 

and John Hovorka. Their book Fire Control Pinciples ^ separates 
the fire control problem from its solution and the principles from 
the actual systems. It deals with "above-water“ fire control and 
assumes the use of electromagnetic radiation to obtain target infor- 
mation. This writer is aware of previous work, some classified, 
in the underwater fire control field, but as far as he has been 
able to determine, this thesis is the first work of its nature. 

It is an extension of the line of attack pursued by Wrigley and 
Hovorka to the surface- to-underwater problem. In some facets the 

two problems are quite similar, but in many there are significant 

(3) 

differences. Concurrently with this research, Messere' * is 
attacking the underwater- to-underwater problem. 

1.3 Outline of Attack 

The problem will be attacked in three broad steps. First 
those quantities not found in the "above-water" problem are defined, 
examined, and illustrated when appropriate. Then the problem is 
formulated following the development in Chapter Two of Wrigley and 
Hovorka. Finally two general methods of solution of the problem 
are proposed, briefly analyzed and compared. 



CHAPTER I I 



qUANimES PECULIAR TO THE UNDERWATER FIRE CONTROL PROBLEM 



2 . 1 I ntroduct ion 

In beginning the investigation of the surface- to-underwater 
fire control problem, one is immediately struck with its points of 
difference with the "above-water" fire control problem. In the 
underwater problemp intelligence as to the target's position is 
received by sound waves which 

(1) travel slowly enough compared to target speeds 

that the above mentioned intelligence must be considered 
as past data; 

(2) travel in a medium where the speed of sound has 
enough variation within the volume of consideration 
that refractive effects cannot be ignored. 

On the other hand, "above-water" electromagnetic radiation can be 
assumed to travel in straight lines for the purposes of the problem, 
and the intelligence of the target’s position it conveys can be 
assumed to be current. 

2.2 Splash Point 

One way of considering the surface- to-underwater problem is 
to determine the position on the surface of the water, called the 
splash point , where a weapon should enter so that it will sink 
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vertically to the future position of the target, or in the case 

of homing weapons to the optimum position relative to the target 

to attack it. Thee geometr i cal features of the splash point are 

illustrated in figure 1. Using this approach the problem is that 

of a surface weapon station delivering a weapon to a given position 

on the surface of the water. The vector R^, splash range , is 

defined as the direction and distance from the present position 

of the weapon station to this given position (the splash point). 

In general, this vector is determined by predicting target motion. 

2„3 The Water-Mass Coordinate Frame 

The concept of target motion brings to mind the following 

question: target motion with respect to what coordinate frame? 

It is of course possible to generalize, but for the purposes of 

this thesis the water-mass frame is used.' This has been chosen 

for several reasons: the weapon station and the target move in it, 

the sound waves which are used to obtain target information travel 

in it, and it can be assumed that a weapon upon entering the water 

at the splash point will fall vertically in the water-mass, i.e., 

that any horizontal component of velocity with respect to the 

water-mass at time of water entry will be essentially brought to 

zero at once. The coordinate frame is arranged in the water-mass 

to yield horizontal and vertical components of vectors. These 

components are designated by a final subscript H or V as 

appropriate. 

- 

See the discussion of coordinate frames in Wrigley and 
Hovorka, pp. 5“7. 
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2.4 Transmission Time 

In the "above-water" prob 1 em, i ntel 1 i gence as to the position 
of the target is usually assumed to be current. Such is not the 
case in the underwater problem. T ransmi ss ion time , denoted by 
the symbol t^, is defined as the interval of time required for 
sound to travel from the target to the weapon station sound detector. 
It is illustrated in figure 2 and quantitatively defined by the 
equation 




where is antecedent range : the distance from the present 

position of the sound detector to the position of the target t^ 

previously' and where v(r) is the speed of sound as a function of 

the path it takes. In practice, of course, an average or standard 

value for the speed of sound is used, for example 1500 meters per 
(4) 

second.' For equal to 3000 meters tj is then equal to 
2 seconds. In this time a 20 knot target will travel over 20 yards, 
a distance that cannot be neglected in the solution of the problem. 

2.5 Line of Sight and Line of Position 

The following quantities are defined in order to formulate 
the problem of determining target motion. Refer to figure 3. 

Ij^g^ - direction of the present 1 i ne of sight : the 

reverse of the direction (as measured at the weapon station) 

'it should be mentioned that R^ is a corrected quantity and 
not raw data. The actual range must be corrected for the refractive 
effects present. 
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of the sound reflected from the target t^ previously. 

l._ " direction of the antecedent line of position: 

LPa ^ 

the direction of the vector from the present position 
of the weapon station (the sound detector) to the ante- 
cedent position of the target (measured t.j. previously). 
The preceding two quantities differ by what is called inaccuracy 
in the line of sight. This inaccuracy is a statistical quantity 
with a mean and a variance. Neither the mean nor the variance are 
known exactly, but they can be estimated in a statistical sense. 
The mean of the inaccuracy is called error and the variations 
about it are uncertainty. Correction to the line of sight , 
symbolized (C)LS, is defined as the negative of the best estimate 
of error in the line of sight. It is determined by consideration 
of the speed of sound versus depth characteristics of the trans- 
mitting medium, assuming that the speed of sound at any one depth 
is constant over the water area involved (see assumption below). 
Now 



thus (C)LS is in the horizontal plane under the following two 



( 2 . 2 ) 



assumpt ions : 

(1) The effects of variation in salinity and pressure 
on the speed of sound are neglected. ' 

Vinsler and Frey^ p. 435. 
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(2) The gradient of water temperature is always vertical. 

2.6 Target Velocity 

Because of the inherent lag in target information caused by 
transmission time, present target velocity cannot be directly 
determined as in the "above-Water" problem. Instead, antecedent 
target velocity , symbolized is the available quantity used. 

It is defined as the direction and speed of motion of the target 
prior to present time. 

2.7 Sinking Time 

Sinking time of the weapon, denoted t^, is illustrated in 
figure 4. It is measured from the time of water entry to the 
time of expected detonation or time of commencing homing as appro- 
priate. Sinking time will be an empirically determined quantity 
dependent upon the weapon and the depth to which it is expected 
to sink. This depth is the vertical component of the future line 
of position, which is defined in section 2.9. 

2.8 Future Range 

Figure 5 illustrates future range Rp. It is the vector from 
the present position of the weapon station to the position of the 
target at future time t^ + t^, where t^ is the time of flight of 
the weapon. Rp is quantitatively defined by the following 
equat ion : 

0 o 



(2.3) 
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where V^(0*) is target acceleration as a function of time. This 
acceleration cannot be predicted except in a statistical sense. 

Let the statistical "best guess" of the average value of V.|.(^) over 
the time interval t.|. + t^ + t^ be symbolized \l^. Then equation 
( 2 - 3 ) becomes 




(2-4) 

2.9 Offset Range and Future Line of Position 

Offset range , Rgpp* is shown in figure 5. It is a vector 
used when homing weapons are employed and places the weapon at an 
optimum target angle and range with respect to the target at time 
tp + tg after firing. The idea of offset range has to do with a 
"weighted" center of lethality. This concept, illustrated in 
figure 6, can be applied to any weapon that does not have to 
strike the target to detonate with lethal results. In the 
case of the proximity fuse the center of lethality is in fact the 
target itself since the velocities of the weapon fragments after 
burst are so much greater than the target's velocity that there is 
no measurable difference in the lethality of a burst equally dis- 
tant ahead of or behind the target. Such is not the case with a 
homing torpedo. After arriving at the end of the offset range 
vector, it homes for the target using its own propulsion, which 
may or may not give it a speed advantage over the target. Clearly 
then the center of lethality for such a weapon is not at the target. 
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but generally ahead of the target. The exact position is deter- 
mined by the characteristics of the weapon and V-j-g- For non-homing 
weapons Rq^p is, of course, zero. 

The future line of position , symbolized LPp, is defined as 
the vector from the present position of the weapon station to 
the time of Rgpp- As is indicated in figure 5: 

LPp = (2-5) 



2.10 Summary: The Determination of Splash Range 

Finally it is seen from figure 5 that since the weapon: 
must sink along LPpy to achieve its objective, the point of water 
entry must be at the tip of LP-u- Thus: 

r H 



R = LP^„ 
s FH 



R 

s 



(2-6) 



R + R 
FH OFFH 



(2-7) 



CHAPTER I I I 



FORMULATION OF THE PROBLEM 



3.1 Introduction: Mi ss-Produc ing Effects 

Chapter II considers some of the aspects of the surface-to- 
underwater fire control problem that differ from the "above-water" 
problem. Now the problem itself is formulated. As noted by 
Wrigley and Hovorka' there are three miss-producing effects: 

(1) Target motion 

(2) Curvature of the projectile's trajectory 

(3) Jump 

In the problem of this thesis, target motion involves the 
travel of the target between its last known position (this position 
being that of the target t.|. prior to firing) and its position at 
the time of expected detonation, or the time of commencing homing, 
as appropriate (this position being that of the target t^ + t^ after 
firing). Unlike the "above-water" problem it also involves a dis- 
placement of the center of lethality from the target's position 
at the time of commencing homing. Since this displacement is 
dependent upon the target's motion (it would be zero if the target 
were stationary) it is included as a miss-producing effect involved 

'wrigley and Hovorka, p. 24. This chapter follows their 
development closely. 
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with target motion. Figure 7 illustrates target motion as a 
miss-producing effect. 

Trajectory curvature is made up, as in the ‘‘above-water" 
problem of the effects due to the forces that act on the projectile 
during its time of flight. It is illustrated in figure 8. 

Similarly, Jump is not altered by the special conditions 
of this thesis and represents the effect of factors that make the 
initial velocity of the projectile differ in direction from that 
of the weapon line, which is def i ned^as the direction in which the 
weapon is aimed. Figure 9 illustrates it. 

3.2 Corrections for Miss-Producing Effects 

If the weapon line were to be oriented along the horizontal 
component of the antecedent line of position the three miss- 
producing effects mentioned above would all cause the projectile 
to miss the splash point. If a perfect set of corrections could 
be developed to apply to the weapon line prior to firing, the 
miss-producing effects would be completely nullified and the 
projectile would enter the water exactly at the splash point as 
desired. Toward this ideal the following corrections for miss- 
producing effects are developed: 

(1) Lead, which compensates for target motion effects. 

(2) Curvature correction, compensating for in-flight 
forces acting on the projectile. 

(3) Jump correction, to compensate for initial 
velocity effects. 



These corrections are illustrated in figure 10. 
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3-3 Prediction Angle and the Prediction Angle Vector 

Since each correction is a directional compensation it is 
therefore an angle. Hence the sum of these component angles is 
an angle and is defined as the prediction angle, measured from the 
horizontal component of the antecedent line of position to the 
weapon line. The vertical component of the antecedent line of 
position does not enter into the definition of prediction angle 
as the problem is formulated as a surface-to-surface problem with 
the projectile being fired to hit the splash point. 

Under certain limitations the prediction angle P can be 
defined as a vector angle.' Consider the cross product of 



1 ! 



LPaH 



with 1 



Wf 



lLPa.H^-^WL’^J.(LPoLH,WL) P 

( 3 - 1 ) 

Now if P is a “small“ angle, i.e., sin P = tan P = P and 
cos P = 1 within limits of allowable tolerances, 
the vector angle P can be defined as 

P - 1 X 1 

' 'kL(LPciH,WLr -^LPaH 

(3-2) 

Since the accuracies of measurable quantities involved in the 
problem of this thesis are not anywhere near those of the 

^Wrigley and Hovorka, p. 19- 
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"above-water" problem, the inaccuracies introduced by employment 
of the "small" angle assumption throughout this thesis will not 
materially effect the accuracy of the solution of the problem. 
Henceforth this "small" angle assumption will be valid, and 
will not produce any restrictions to the problem. 

3.4 Lead 

Lead is defined as the angle between the horizontal component 

of the antecedent line of position and the horizontal component 

of the future line of position. Unlike the "above-water" problem, 

the future line of position is the direction from the present 

position of the weapon station to the point below the splash 

point defined by the tip of the offset range vector. Thus the 

future line of position can be thought of as pointing toward a 

pseudo-target . Note that 1 ,dcli identical to R /R . Figures 

Lrrrl S S 

5 and 11 assist in understanding the following derivation: 

Under the ’’small" angle assumption 



L = 1 x1 

^ -^LP<xH -^LPFH 



1 



LPouH 



R, 



R, 



(3-5) 



Now 



(1) Rg Rp^ + Rqppi^ 

(2) P r \/ t 

where average projectile velocity 



(2-7) 



(3-6) 



-14- 



The horizontal component of eq. (2-4) gives 
(3) 



(3-7) 



^ ^ ^LPgH^^OPPH 



Vp(av)tf 






(3-8) 



Now 



/ V X ^ V 

(W \ ~ — LPq,H ^Tc^H 

V''\paH^T P 

rx^kU 



(3-9) 



where angular velocity of the horizontal component 

of the antecedent line of position that is due to target motion. 

^LPaH "" ^^LPaH^T ^\paH^WS (3-10) 

where is the angular velocity of the horizontal component of 

the antecedent line of position and is the angular 

velocity of the horizontal component of the antecedent line of 

position that is due to weapon station motion. 

1 



Wrigley and Hovorka have shown that 

/—V I 

(W \ =--±kP^±i_W 

\ ''LPclH7ws P 



OlH 



(3-11) 



Wrigley and Hovorka, p. 32, figure 2-6, 
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where weapon station velocity, 

Thus, by combining and rearranging 



I _ RaH 

pCw) 






V, 



+ ^\PaH V 

v/L Y c 



'pCav) 






(tr*tf4sf , Ilpcm^Ro 

. .j 






zv, 



p(OLV) 






^P(olV) 



( 3 - 12 ) 



it must be noted that !,» u cannot be measured, but 1.. ,, 

LPaH LSoH 

can. Fortunately, (C)LS is in the vertical plane under the 

assumptions made in the last chapter (see the section on line of 

sight and line of position) so l,n „ and l,e u identical. 

Similarly, while W, ^ cannot be measured, it is the same as 

LPaH 

W ^ the angular velocity of the horizontal component of the 
present line of sight, which can be instrumented. 

3-5 Curvature Correction 

Curvature correction is illustrated in figure 12 and 
is defined as a rotation from the horizontal component of the 
future line of position to the projectile line. Thus, under 
the "small" angle assumption the curvature correction vector is 
defined as follows: 



^ ~ 'iPFH ^ ’pL 



(3-13) 



where C is the curvature correction 
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1 . = Vr. /Vr. > the unit vector along the projectile line 

PL Po Po ^ 

Vp^ is the projectile velocity at the time of firing 
Splash range is given by 



o O (3-14) 

where Vp is the projectile acceleration 
Splash range is also given by 

% " \pFH *^s " \pFH''p(av)’^f (3-15) 

Substituting equations (3-14) and (3-15) into equation (3-13) 
gives 




^p(O.V^ 



( 3 - 16 ) 

For the purposes of this section of the thesis the 
water-mass is assumed to be effectively inertial space and thus 

" ^P (3-17) 

where fp is the sum of the principal specific forces (forces 
per unit mass) acting on the projectile during its flight. 

These principal forces are: 

(1) Gravity, G, the gravitational field intensity due 
to the mass attraction of the earth for the projectile. 
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(2) Aerodynamic drag per unit mass; for Mach numbers 
encountered in the problem of this thesis, this is 

” '^(AM)P '^(AM)P’ empafical coefficient 

and is the projectile velocity in the air mass. 

(3) Lift per unit mass, f^j,^. 

(^) Thrust per unit mass, when the projectile 

i s sel f-propel 1 ed . 



Thus , 



'^P ^ ~ ^0 '^(AM)P '^(AM)P ^(li) ^(th) 



(3-18) 

Integration of this equation involves line integrals, as the 
vectors are generally changing over the time interval of the 
integration. Howeverj using average values of the vectors over 
the t ime interval : 



rV- f ^ 

{{y/r)dvdt=i 



"ft 



G ko(av^ 



^ ^(nXo.v') 



Substitution of equation (3-19) into equation (3“l6) gives; 



■f f X 1 



F 

PL '(thXov) -^PL 



( 3 - 20 ) 



It should be noted that no drift term is present since ASW 
projectiles have no significant spin. 
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3-6 Jump Correction 

Figure 13 illustrates jump correction, which is the angular 
correction needed to compensate for the non-parallelism of the 
weapon line and the initial projectile velocity. It is measured 
from the projectile line to the weapon: line. The jump correction 
vector is the negative of the jump vector, J, and under the 
"small" angle assumption is defined below: 

- J = »pL X l^L (3-21) 



For the problem of this thesis, - J is made up entirely of 
weapon station velocity jump correction. It is required because 
a projectile leaving the weapon station with a velocity '^[(ws)P]o 
relative to the weapon station has the total velocity: 



'^Po '^[ (WS)P]o '^WSo 



( 3 - 22 ) 



where the weapon station velocity at the time of firing. 

Windage jump' is not present as ASW weapons have no signifi- 
cant spin. 

is along 1 



'pL ' ''PO^PO ''[(US)Plo 



WL 





(3-23) 



3.7 Summary 

Prediction angle, P, is defined as a vector rotation from 
the horizontal component of the antecedent line of position to 

'wrigley and Hovorka, p. 39. 
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the weapon line. Under the "small" angle assumption the following 
equation is valid: 



^ ’iPaH ^ ’wL 



(3-2) 



Thus, prediction is given in terms of lead, curvature 
correction, and jump correction as: 

P = L + C - J (3-24) 

Upon substitution of the previously derived equations (3-12), 
(3-20), and (3-23) > prediction angle is seen to be made up of the 
fol lowing terms: 












Vp(av) tf 



LPaH 



V, 



PCotv) 






4- 



Horizontal component of line of Weapon station 
position angular velocity velocity 



I I 1 

-*-LPaH I 4 J f 



>^Lead 



2V 



P(Q.n.> 






V t 

T(a^) '-f J 



Horizontal target acceleration Horizontal offset 



L 



r 



zy 



P(aN») 



^ ^ 1 PL ~ P(o.'»)Vam^ P(iav) ^ 1 Pl 



Gravity drop 



Aerodynamic drag 



-f- 



f^(li)(ov) f*( 






L] 



4- 



Curvature 

Correction 



Lift and thrust 



4 “ 



VvJSo'^ 

Vfo 

Veloc i ty j ump 



Jump 

Correct ion 

(3-25) 



The entire development is summarized in derivation summary 3“ 1 . 



DERIVATION SUMMARY 3“ 1 





P = L + C - J 


(3-24) 


By 


def in i t ion : 






1” 5 T X T 

LPaH LPFH 


(3-5) 


But 


^LPFH " 




and 


" ^FH ^OFFH 


(2-7) 




^F = "a ^ ^Ta ^^T ^ 'f ^ 





+ \ (2-4) 

Taking the horizontal component of R^ 

^FH ^ ^aH ^TaH ^^T 2 '^TH 

(3-7) 



For the projectile to hit the splash point 

'^p(av)^f (3-6) 

So by combining and substituting 




luft-W s) I LPoiH^ 
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Now 
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LPaH 



^'*^LPaH-^WS 
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By definition: 
C = 



'lpfh X Vl 
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Now 
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Thus 



7- 

^ V tc 

''PCqv) 



PL 



( 3 - 16 ) 

Now 

« — 

'^P ~ (3-17) 

^ ■ *^D ''(AM)P '^(AM)P "^(1 i) ^(th) (3-18) 

Using average values 
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By ( 


def i n i t i on : 






- J = "i’pL ^ ”^WL 
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But 
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Final ly, from equations (3"24) , (3"12), (3”20) and (3”23) 



R.h (tAA) y , 



V, 



P<0'») 



t 



LPaH W 






IlPolH ^ RopPH 



zv. 



PCov) 












2V, 



rc«) 



^ ^ ^ PL J 







\WL 




(3-25) 



CHAPTER IV 



SOLUTION OF THE PROBLEM 



4. ] I ntroduct ion 

The preceeding two chapters have defined certain quantities, 
given a kinematic and geometric look at the problem, formulated 
the problem, and finally derived an expression for prediction 
angle. If all the quantities involved in this expression could 
be measured or computed (or even statistically "best-guessed” 
at) the problem would be solved. The solution, then, revolves 
around these quantities. 

Many of the quantities involved are obtained in the same 
manner as in the "above-water" fire control problem. However, 
attention will be focused on those quantities that are peculiar 
to the underwater problem. All of the peculiarities occur in 
the lead part of the prediction angle and come about primarily 
as a result of using relatively slowly traveling sound waves to 
determine the position of the target. In addition, the quantity 
t^ is present as a consequence of the target's being below the 
surface of the water. 

4.2 Information Rate 

The use of sound causes an appreciable transmission time, 
as discussed in chapter II, a typical value being two seconds. 
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Furthermore the relatively low velocity of sound, about 1500 meters 
per second, means that sound pulses must be sent out at well 
separated intervals of time rather than semi-cont inuous ly as in 
radar. For example in radar 800 pulses per second is a typical 
rate^^^ giving a maximum unambiguous range of slightly less than 
200 kilometers, while for a maximum unambiguous range of 7500 
meters using sound waves in water the interval between pulses 
must be ten seconds. Under the conditions of the example the 
best information rate possible for the underwater target is one 
piece of information every ten seconds. Figure 14 illustrates 
the concept of this section. 

4.3 Transformation of Target Information 

Another complication comes about from the fact that it is 
not practical to stabilize the sound transducer against the 
angular motion of the weapon station. This means that informa- 
tion as to the location of the target is available only referred 
to the structure of the weapon station. A coordinate transforma- 
tion must be used to transform target information into a frame 
that can be instrumented and related to the water-mass frame, in 
which the computation of the prediction angle is to be carried 
out. An example of this process is described in Appendix A. 

4.4 The Trianqulat ion Method 

If successive positions of the target are used to determine 
lead, the following expression can be used: 
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\ _ luFbLH^RopPM 

L^ ■*“ 



Vc 



'P<av7 tf 

Horizontal target motion 
Horizontal Target Acceleration 






P<av> 

Horizontal offset 



(3-8) 



The quantity is obtained by coordinat izing target positions in 

the water-mass frame (see fig. 15)- By subtracting successive target 
positions its change in position relative to the weapon station 






[ (WS)T] 



can be calculated. Since the weapon station velocity 



and the :ime between target positions are known, the change in wea- 
pon station pos i t ion , , can be calculated and combined with 

the previously computed change in target position relative to 
the weapon station. This gives the change in target position in the 
water-mass hK- The horizontal component of this change divided by 

r 

the time between target positions, /^t, is and the vertical com- 
ponent of the change divided by the same time interval is Thus, 

V.^. is computed. This quantity, together with the type of weapon used, 
T a 



determines Rgpp and hence 



V” is not directly available from 



target information. It might be computed by dividing the change be- 
tween successive values of V_ by the time interval between them and 
assuming a constant acceleration over this interval. It also could 
be set in to the computer by hand. Transmission time, t.^., is empirically 



determined. The quantity 'S, of course, available as 
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and all the other quantities in equation (3-8) are determined as 
in the "above-water" problem. Thus, lead can be computed. Curva- 
ture correction and jump correction deal with "above-water" 
effects and are computed as in the "above-water"problem. The 
method of solution of this section is called the triangulation 
method in this thesis. 

4.5 The Controlled Line Method 

In this section the original expression for lead, equation (3”12), 
is used. The quantity must be measured in order to get 

As noted by Wrigley and Hovorka', the line on the fire control 
equipment that tracks the line of sight is called the tracking line, 
or indicated line of sight, and angular velocity of this line is 
the best information available as to the quantity Since it 

is not feasible to stabilize the sound transducer, the method of 
this section could use either of two approaches: 

(1) Use a stabilized dummy tracking line and an 
unstabilized tracking line. 

(2) Transform the target information prior to displaying it. 
Using the first approach, the components of dummy tracking 

line orientation, expressed in the reference coordinate system of 

Appendix A, would be transformed by the inverse of the matrix of 

equation (A-8) into components in the weapon station structure. 

These components would position a tracking line on displays of 

2 

unstabilized target information. A lack of coincidence between 



'wrigley and Hovorka, p. 53 
2 

Here coincidence is used in the sense of servo null: reducing 

the difference to the smallest value commensurate with the servo 
configuration, not necessarily to zero. 
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this tracking line and the displayed line of sight would give rise 
to a tracking inaccuracy« Thus, a tracking inaccuracy correction 
would have to be generated. This correction is denoted (C)TL 
and is called correction to the tracking line. Illustrated in 
figure 16, it is the angle from the tracking line to the line 
of sight. In the process of solving the problem it would be 
nulled' and thus would be a “small" angle. The vector represen- 
tation of this angle, (C)TL, would be coord i nat i zed in the weapon 
station structure. As shall be seen shortly, it would have to be 
transformed into water-mass coordinates using the matrix of 
equation (4-8) before it would be useful in solving the problem. 

Using the second approach, the target data would be trans- 
formed into water-mass coordinates prior to display. The display 
equipment, although physically a part of the weapon station structure, 
would model the water-mass frame. The orientation of the dummy 
tracking line would directly position the tracking line of 
the displays and (C)TL would be coord i nat i zed in the water-mass 
f rame. 

Since the water-mass frame can be related to inertial space, 
an attractive means of stabilizing the dummy tracking line is pre- 
sented by the space integrator. This device is discussed in detail 

(G ^ 

by Draper, Wrigley, and Grohe' . In essence, the space integrator 
performs two functions: (1) With no command input it holds the 

direction of a controlled line constant in inertial space in the 

'here coincidence is used in the sense of a servo null: 
reducing the difference to the smallest value commensurate with 
the servo configuration, not necessarily to zero. 
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face of interferences; (2) With an input command it gives the con- 
trolled line an angular velocity with respect to inertial space 
that is proportional to the input. So if some input to this 
device can be found that will enable its controlled line to track 
the line of sight, the controlled line will be the dummy tracking 
line and the input command to the device will be proportional to 
the angular velocity with respect to inertial space of the line 
of sight. The horizontal component of this quantity is then the 
desired quantity There are two questions to be answered: 

(1) Can the water^rnass frame be considered effectively inertial? 

It can, as the small errors introduced by this assumption are 
insignificant compared with uncertainties in the line of sight, 
for example: The line of sight is assumed for the sake of general- 

ity, to have an uncertainty of 10 mi 1 1 i radians, which is a realistic 
figure. The interval between successive bits of target information 
is 10 seconds. Then the uncertainty in is 2 mr/sec. Assuming 

that the rotation of the water-mass with respect to the earth is 
of lesser magnitude than earth rate, i.e., the problem is not 
being worked near a whirlpool, then the difference between angular 
velocity in the water-mass and inertial angular velocity is no 
greater than earth rate, which is 0.0729 mr/sec. As this differ- 
ence is much less than 2 mr/sec, the assertion is correct. The 
second question to be answered is: (2) What is e.. suitable. 



command input to the space integrator? Since the space integrator 
operates with respect to inertial space, the input to it must be 
expressed in inertial space. Correction to the tracking line. 
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expressed in the water-mass frame, which has been shown to be 

effectively inertial, is a suitable input when properly modified. 

It will drive the controlled line so that for a constant W, „ , the 

tracking line will be coincident with the line of sight each time 

the line of sight is displayed. As an example, modify the signal 

proportional to (C)TL by multiplying it by R and passing the 

0 

product through a parallel integrator and bypass. Add the output 

2 

of the parallel elements, divide by R , and use the resulting sig- 

0 

na) as the input command to the space integrator (see figure 17 ). 
If the system starts with an initial (C)TL, the integrator will 
build up a charge while (C)TL is being nulled. For a constant 
this integrator charge will give the correct signal to drive 
the space integrator at equal to with (C)TL nulled. 

The favorable answers to the two questions above indicate 
that the problem can be instrumented in this manner. With W 



LPaH 



obtainable from W, ^ j the quantities t_, t , 1,. ,,, and V_ 
LSo’ ^ T s’ LPaH T 



would enter the problem solution as in section 4.4. With the 

exception of R^rru other quantities in equation (3"25) would 

be determined as in the "above-water” problem. However, to get 

V_ must be known. It would be determined essentially as 
OFFH Ta 

in section 4„4 This indirect requirement for is a disadvantage 
of this method. 

In "above-water" systems the instrumentation of the tracking 



line as the controlled line of a space integrator is common. 
Electromagnetic radiation is used and (C)TL can be considered as 
always available. This is not so in the underwater problem. 
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(C)TL is available only each time line of sight information appears 
on the displays. This occurs only when a sound pulse is received. 
Thus, this method of solving the problem involves sampled data. 

The charge on the integrator discussed earlier in this section 
acts as a memory and drives the space integrator between bits of 
data, but a measure of the coincidence of the tracking line with 
the line of sight is available only with the receipt of a sound 
pulse from the target. 

The method of solution employing the instrumentation 
described in this section is called the controlled line method 
in this thesis. 

4,6 A Simple Analysis of the Preceding Two Methods 

How well do the preceding two methods of solution do the 
job of instrumenting the perfect set of corrections mentioned in 
section 3-2? This analysis will compare the propagation of un- 
certainties in target position through the two methods, since the 
methods differ only in their computation of lead. Specifically, 
the following two expressions that give target motion (less accel- 
eration) must be compared: 

(1) From eq. (3"8) 



J-LPaH 
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(2) From eq. (3“ 12) 



R 






V. 



(tM) ^ 

^LPaH^ 






P<av) 



tf 



V 



P(av) 



t, 



The other terms in equation (3“8) also appear in exactly the 



t_ + t, + t 
T f s 



same form in equation (3“12). Further, since 

P(av)^f 

is common to both expressions, the comparison is reduced to 



that between l^p^^ X and X 

'^LPaH = ('^LPaH^ ^'^LPaH^WS 



( 3 - 10 ) 



So 

But 



So 



^'^LPaH^WS 



*LPaH ^ ^WSo 



’^an'^LPaH 'iPaH ^ '^WSo ’^aH^'^LPaH^T 



^'^LPaH^T 



’iPafi ^ ^TaH 



'^aH^'^LPaH^T 'bPaH ^ '^TaH 



(3-11) 

(4-1) 

(3-3) 

(4-2) 



Finally, from equations (4-1) and (4-2) 



'iPaH ^ '^TaH '^aH'^LPaH ’iPaH ^ ^VISo 



(4-3) 



In the controlled line method, is dependent upon input bits 

of information, not a continuous flow as in “above-water" systems. 

However, if is assumed constant, the signal modification 

scheme outlined in the preceding section contains a constant 
2 ~ 

quantity: kR^ '^LSo k is a constant of proportionality. 

This was first observed by Kenneth Wallace about ten years ago. 



Since the quantity R can be generated by the "memory" of the 

a 
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computer between input bits of information, can be kept up 

to date and the W, „ „ determined by the controlled line method 
will correspond to the latest antecedent line of position. Thus, 
equation (^-3) is valid for a comparison of the two methods, and 
indeed shows them to be effectively the same. Thus, the uncertainty 
in lead caused by uncertainty in target information will be the 
same for both methods. A quantitative measure of this uncertainty 
is given by the following illustrative example: 

Assume a non-homing weapon is used and target acceleration 
is zero. 

Let R = 2000 yards and R u “ 1992 yards 

a an 

Then t^ = 1 second as the speed of sound is assumed here 

as 2000 yards per second. 

Let V / X = 1600 feet per second 

p(av) 

t^ = 4 seconds 

t =5 seconds 
s 

Let uncertainty in the horizontal component of the line of sight, 

(U)LS M = 5 m i 1 1 i rad i ans . 

On 

Let maximum unambiguous range be 5000 yards, giving an information 
rate of 1/5 bit per second. 



Then 




= 2 mr/sec - (U)W^p^^ 


Let 




20 yards 




<")''uso 


LPaH ' ' 




“wso UaH " 




w 

LPaH 


10 mr/sec 
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Now under the assumptions 



tp Hp.h tp 






(4-4) 
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Since is not «l , conventional methods of error analysis 
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not appf*i*cable and a straightforward computation is used. 
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I = 1 

Defining lead uncertainty, (U)L, as follows: 





(^- 5 ) 



(U)L = j (.162 - .119) 1^ 



.022 1 



V 



(4-6) 



Both methods are classified as TT systems under the scheme 
of Wrigley and Hovorka'. Although it would be possible to con- 
figure them otherwise, there is no reason for so doing. Indeed, 



in general, better than a WW system, assuming independent operation 
of axes. 

4.6 Summary 

This chapter considers the solution of the problem. The 
sampled data nature of target information is pointed out, as is 
the required transformation of this information. Two methods for 
solving the problem, the t r i angu 1 at ion method and the controlled 
line method, are described and the methods analyzed. This 
analysis finds them to be effectively the same. Equation (4-3) 
relates the quantities that differ in the two methods. 

It would seem then, that there is little to choose between 
the two methods with regard to accuracy considerations. The 



Markey^^^ has concluded that the three dimensional TT system is. 



Wrigley and Hovorka, p. 57- 
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control led line method does have the disadvantage of requiring 

an additional computation of V , needed to determine when 

1 3 OFF 

homing weapons are used. 

An illustrative example of the propagation of uncertainties 
yields an uncertainty in lead of 22 milliradians under the 
assumed conditions. 



CHAPTER V 



SUMMARY AND CONCLUSIONS 



5.1 Problems Caused by the Target *s Being Underwater 

The use of sound to obtain target information is the big 
difference between the “above-water*' and the surface- to-underwater 
fire control problem. This thesis has shown in the example of 
section Z.k that a finite transmission time of sound from the 
target to the weapon station requires the use of past data to 
solve the problem. Concomitant with this is the concept of a 
slow information rate, leading to having some data available 
only as sampled data with a sampling rate which cannot be 
increased beyond a fixed limit set by unambiguous range con- 
s i derat ions , 

The fact that water temperature variations commonly en- 
countered at sea cause the speed of sound in water to vary sig- 
nificantly, requires corrections for refractive effects. The 
footnote on page 5 and section 2.5 mention two of these corrections. 
Fortunately, the angular correction to the line of sight lies 
essentially in the vertical plane, a fact which is of assistance 
in the instrumentation of the problem. 

Although not associated with the use of sound, the concept 
of an offset center of lethality, caused by the fact that the 
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speeds of homing weapons are of the same order as target speeds, 
is used in formulating the problem of this thesis. This concept 
leads to the development of an offset range vector (section 2,9) 
Final ly^ sinking time, t^, not found in the “above-water” 
problem, was used to account for the interval of time from water 
entry of the weapon at the splash point to the arrival of the 
weapon at the point of detonation or the optimum attack position 
as appropriate. 

5«2 Methods for Solving the Surface- to-Underwater Fire Control 
Problem 

Two methods for solving the problem are discussed in this 

thesis: the triangulation method and the controlled line method 

A brief analysis indicates that both methods compute prediction 

angle in essentially the same manner, using the same inputs to 

develop essentially the same quantities by slightly different 

means. Thus, the two methods are shown to be effectively the 

same (equation 4-3) so that inaccuracies in input information 

propagate in the same manner in both methods, leaving no 

choice between them on that account. In section 4.6 it is noted 

that the controlled line method has the disadvantage of an 

additional required computation for V_ when homing weapons 

T a 

are used. 

5.3 Suggestions for Further Work 

Since this thesis is an initial attack on the surface-to- 
underwater fire control problem from the viewpoint expressed 
by Wrigley and Hovorka, more work remains to be done. An 



investigation of the sampled data aspects of the 
one area. 

A frequency diversity scheme might be used 
information rate. This would involve "hardware" 



problem is 

to increase the 
modifications. 



The feasibility of the idea might be investigated. 
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Weapon Station Splash 

Delivers Homing Weapon Point 





Figure 1. Geometrical Features of the Splash Point 
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Target When Sound Reflected, 
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The Relationship Between Line of Sight and Line of 
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Figure 3- 
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at Splash Point (time t + t ) 

(time t) ^ 

a) The Homing Case 




at Splash Point (time t + t ) 

(time t) ^ 

b) The Non-Homing Case 



Figure 4. Sinking Time 
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Figure 5. The Relationship of Future Range and Offset Range 
to Splash Range 
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Figure 6. Center of Lethality 
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Figure 8. Curvature of the Trajectory as a Miss-Producing Effect 
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Figure 10. Geometrical Features of the Complete Problem 



- 47 - 



Splash Range 
R 



s P(3v) f 

for Projectile 
Impact at 
Splash Point 



Lead Angle 
Vector 



Horizontal Component 
of Future Line 
of Position 




ash 
Point 

Angular Velocity 
of Horizontal Component 
^of Antecedent 
Line of Position 



Hor i zonta 1 
Target Motion 
Composed of 
Horizontal 
Jarget Travel 
During Interval 

(t^ + + tg) plus 

Horizontal Offset 



' Point Directly 
\ Above Antecedent 
Posi t ion of the 
w Target 



Weapon Station 
Position at 
Time of Firing 



zonta 1 Component 
Antecedent Line 
of Position 



Horizontal Component 
of Antecedent Range 



R 



aH 



Figure 11. Geometrical Features of Lead Angle 
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Figure 12. Geometrical Features of Curvature Correction 
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Figure 13. Geometrical Features of Jump Correction 



- 49 - 
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Figure 14. 



Information Rates in Electromagnetic and Sound 
Rad iat i on 



- 50 - 




V 

V 

V 



ta 

TaH 

TaV 






Ta 



Figure 15- Determination of V 
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APPENDIX A 

AN EXAMPLE OF TARGET INFORMATION TRANSFORMATION 



Consider the situation where the roll angle, 9, the pitch 
angle cf, and the heading, of the weapon station are obtained 
as fol lows. ' The origin of the reference coordinate system is 
fixed in the weapon station and it has x^, y^, and axes 
oriented north, east, and down respectively. The weapon station 
structure is initially aligned with the reference system. Con- 
sider figure A- 1 . A heading of 4^ aligns the structure with the 
axis system Xj, y^ and z^. Note that since heading is measured 
about the z^ axis that the z^ and z^ axes are coincident. 

Pitch angle is measured relative to the Xj , y^, z^ 
system and is developed by a rotation of the structure about the 
yj axis into the x^s ^2^ ^2 system. Note that the y^ and 

y^ axes are coincident. 

Finally, from its alignment with the x^ , ^2^ ^2 system, 

the weapon station structure rolls about the x^ axis into align- 
ment with the x^, y^, z^ axis system. Here the x^ and x^ axes 
are coincident. 

^ The sequence used here is employed because a gyrocompass 
can measure 4^ and because most weapon stations roll to much 
greater angles than they pitch. This writer has experienced 
roll angles in excess of 50 degrees in a destroyer. 
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The weapon station receives target information that can 
be represented by a vector R. Its coordinates are known relative 
to the weapon station structure, i.e., in the x^, y^, axis 
system. The transformation to the reference coordinate system of 
Xq, y^, Zq is conceptually done in three steps. Let 
symbolize a column matrix representing the vector R coord i nat ized 



in the x , y , z axis system, 
n ' n n ' 

Then [R ]2 = [A][R]^ 

[R], = [B][R]2 = [B][A][R]3 



[R]q = fC][R], 


= [C][B][A][R] 
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seen fro.m figure 
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0 cosG 
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sin^ 
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cos<(> 
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cos T 
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C 05 V 


0 
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By matrix multiplication 




COS^COSH^ 

-COSCjiSinH^ 

-S>n4 



Si•^T 

COS? 

0 



[C]tB]lA]- 



Sin <l> Cosy 
-sin^sinT 

C0S(^ 



(Ar7) 



COS^COS^P C05^SlnH^^Sin©5/n4>cos4^ SiH© S‘n4^+CCS^^'V^^cos4*; 

/ 

-^COS(pSln^ C0S((>C0Sf’^ S^ne$;n4> S'OM" S*^6cos4^-cos05ir^45.'n'f / 

-Sin^co$<(t Cos^cosp ^A-8) 

Employing the transformation matrix of equation (A-8) , 
target information is transformed from the axis system of the 
weapon station structure into the reference coordinate system. 

The reference coordinate system of this example is the same as 
the water-mass frame for coord i nat i z i ng position vectors. Since 
the weapon station is moving relative to the water-mass, this is 
not so for velocity vectors. However, the target information 
available is all position data, so the transformation of equation 
(A-8) gives this information coord i nat i zed in the water-mass 



frame. 



■T ■ 
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« • I 









I 



1 
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From initial alignment in the x^, y^, frame, the weapon 
station structure is carried into the , y^, Zj frame by 
heading, *4^, thence into the x.^, ^2 pitch, and 

finally into the x^, y^, z^ frame by roll, 9. 







Figure A- 1 Axis Systems in Weapon Station Structure Angular 
Motion 
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